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Abstract. The purpose of this work is twofold. The first is to survey some parametriza- 
tions of positive matrices which have found applications in quantum information 
theory. The second is to provide some more applications of a parametrization of 
quantum states and channels introduced by T. Constantinescu and the last author, 
and thereby to provide further evidence of the utility of this parametrization. This 
work is dedicated to the memory of our colleague and teacher, the late Professor T. 
Constantinescu. 



1. Introduction 

Positive matrices play a vital role in quantum mechanics and its applications (in 
particular, quantum information processing). Indeed the two basic ingredients in the 
theory of quantum information, viz., quantum states and quantum channels involve 
positive matrices. See, for instance, [H] EI]. Thus, a study of parametrizations of 
positive matrices seems very much warranted. In particular, the very useful Bloch 
sphere picture, [HI EI], for the quantum state of a qubit has prompted several attempts 
at the extension of this picture to higher dimensions. In the process, several groups 
of researchers have looked into the question of finding tractable characterizations of 
positive matrices, which could lead to useful parametrizations of positive matrices, 

This paper is organized as follows. In the next section we set up basic notation 
and also point out some sources for positive matrices in quantum mechanics and its 
applications. The third section introduces six (perhaps well-known) characterizations 
of positive matrices, and reviews some putative parametrizations of states of qudits. 
In the next section, we review a parametrization proposed in [5], reiterating its utility. 
The final section offers two more applications of the parametrization in ||. The first 
concerns Toeplitz states, i.e., density matrices which are also Toeplitz. The second 
investigates constraints imposed on relaxation rates of an open quantum system by the 
requirement of complete positivity. 



2. Sources of Positive Matrices in Quantum Theory 

Let us recall that a matrix is positive semidefinite (positive, for short) if z*Pz > for 
all z G C n . One can easily extend this definition to infinite positive matrices. In effect 
such a matrix is what is called a positive kernel, [H], viz., a map K : iVo x N Q — > C, 
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where Nq is the set of non- negative integers, with the property that for each n > 0, 
and each choice pi, . . . ,p n in N and each choice Zi, . . . , z n of elements of C we have 

n 

^ K(p i ,p j )z i z j > 
».i=i 

Positive matrices intervene in at least two of the basic ingredients of quantum me- 
chanics and quantum information theory viz., quantum states and quantum channels. 
There are, of course, more sources for positive matrices, but, for reasons of brevity, we 
will confine ourselves to discussing states and channels. 

The state of a d- dimensional quantum system is described by a d x d positive density 
matrix of trace 1, that is, a positive element of trace 1 in the algebra Aid of complex 
d x d matrices. States described by rank one density matrices are called pure states. 

A quantum channel is a completely positive map $ : A — > C(TC) from a C*-algebra 
A into the set C(Tl) of all bounded linear operators on the Hilbert space Tt (in the 
situations most frequently met in quantum information processing, A = Aid and 
C(TC) = Ai d >). By the Stinespring theorem, [19J, Theorem 4.1, such a map is the 
compression of a *-homomorphism. For A = Aid, there is a somewhat more explicit 
representation, given in jlj (see also [IB])- Thus, $ : Aid — > £>(T~t) is completely positive 
if and only if the matrix 

(2-1) 3 = S* = mE kjj )} d k> . =1 

is positive, where E^j, k, j = 1, . . . , d, are the standard matrix units of Aid- Each Ekj 
is a d x d matrix consisting of 1 in the (k,j)th entry and zeros elsewhere. 

Remark 1. Usually one requires a quantum channel to satisfy two additional require- 
ments: i) $ be trace preserving, and/or ii) $ be unital. 

A Kraus operator representation of a completely positive map is a (non-unique) 
choice of operators V{ such that one can express the effect of $ via 

r 

$( P ) = Y,v iP v; 

i=l 

Usually only the non-zero Vi are taken into account in the above equation (though 
sometimes it is convenient to ignore this convention). 

Then $ is trace-preserving iff Yli=i = while $ is unital iff Y^i=i = ^d. 
These properties can also be verified (without any reference to Kraus representations) 
by computing the partial traces of 5$ viewed as an unnormalized state (see 

All choices of Kraus operator representations for $ come from square- roots of S$, 
i.e., matrices T such that S$> = TT*. One then obtains the Vi from the ith column 
of T by reversing the vec operation, [TlH I26j. Recall that the vec operator associates 
to a d x e matrix, V, a vector in C de obtained by stacking the columns of V. It is 
precisely because of lack of uniqueness in the square roots of S$ that the Kraus operator 
representation of $ is non-unique. 
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We should point out that some of the definitions for quantum channel notions used 
in jSJ, though equivalent to the standard ones (i.e., the ones used here), are different. 

3. Characterizations of Positive Matrices 

All positive matrices are Hermitian (unlike the real case, the definition of a posi- 
tive matrix automatically forces Hermiticity). There are several characterizations of 
positive matrices as a subclass of Hermitian matrices. Some of these yield useful 
parametrizations of positive matrices. 

The following theorem, which for the most part is standard textbook material (see, 
for instance, the classic ^1]), reviews some of these characterizations. 

Theorem 2. Let P be a Hermitian matrix. Then the following are equivalent: 

• PI P is positive. 

• P2 All the eigenvalues of P are non-negative. 

• P3 There is an upper-triangular matrix T such that P = T*T ( Cholesky de- 
composition) 

• P4 All principal minors of P are non-negative. 

• P5 Letp{t) = t n + Y^i=o (— l) 4 M ra_ * be the characteristic polynomial of P. Then 
bi > 0, for all i. 

• P6 There is a Hermitian matrix H such that P = H 2 . 

Remark 3. i) P3 is normally mentioned only for positive definite matrices in the 
bulk of the literature. However, a limiting argument shows that it is valid for positive 
semidefinite matrices as well, ii) That P5 is equivalent to P2 is just a consequence of 
Descartes' rule of signs. Hi) P4 should be folklore. Quite surprisingly, we were unable 
to find any source where P4 is stated explicitly (even in a venerable text such as [TT]). 
Since a similar statement for positive definite matrices (viz., positive definiteness is 
equivalent to the positivity of the leading principal minors) is well documented and 
we have seen this statement occassionally incorrectly applied to positive semidefinite 
matrices, we will include a brief proof here. Clearly if P is positive, all principal sub- 
matrices of P are positive, and hence all principal minors are non-negative. Conversely 
suppose all principal minors of P are non-negative. Since the coefficients bi of the char- 
acteristic polynomial of any matrix are just the sum of all the i x i principal minors of 
P, it follows that bi > 0. Hence P is positive. 

Whilst the above conditions are equivalent to positivity, they typically do not lead 
to useful parametrizations of positive matrices. For instance, parametrizing P by its 
eigenvalues only describes the U (n) orbit to which P belongs. For the same reason one 
cannot parametrize P by the coefficients bi of the characteristic polynomial p(t) . 

However, one can turn these characterizations into potential parametrizations. To 
illustrate this consider the problem of parametrizing quantum states in dimension d, 
i.e., d x d positive matrices with unit trace. The standard starting point is to represent 
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a state p via 

d 2 -i 

(3.1) p=-(I d +J2fc\) 

i=l 

Here $ G R and the A« form an orthogonal basis for the space of traceless Hermitian 
matrices, specifically Tr(AiAj) = 2^. One typical choice is the so-called generalized 
Gell-Mann matrices, |Tot l2"7j . This basis is obtained from the matrices Ej.j,k,j = 
1, . . . , d (E k j = eke*) via the following construction: 

flj = E ktj + E jik , k < j, 
ft = ~ ( E j,k ~ Ekj) , k>], 



hf = I d , ht = ht 1 ®^ I <k <d, hi = J-^j-^(ht 1 ®(l-d)). 

These matrices, fkjjhf,hf.,h^ together form one choice of the {Xi,Id} basis for the 
space of d x d Hermitian matrices. When d = 2 this is precisely the Pauli matrix basis. 
When d = 3 one gets the usual Gell-Mann matrices. 

With Equation ()3.1|) as the starting point one can restrict the vector /3 — (fli, . . . , b^^] 
.R d2_1 to satisfy any of the characterizations PI - P6. In principle, this provides a 
bijection from a subset of R d2 ~ l , say Dp, to the space of d x d density matrices. This is 
precisely what is proposed simultaneously in 0] for the characterization P5. How- 
ever, now by conservation of difficulty, the burden of the analysis of quantum states in 
dimension d > 2 is shifted to obtaining a concrete analysis of the subset Dp. In par- 
ticular, these do not lead to easily computed parametrization of quantum states (cf., 
the conclusions section of |T3] ) . Interestingly enough each of these characterizations 
leads precisely to the Bloch sphere picture when d = 2, as we encourage the reader 
to verify. However, this approach has some utility in higher dimensions as well. For 
instance, depending on which characterization one uses, it is at least possible to be 
more precise about the set of pure states (i.e, rank one states). We shall explain this 
via the characterization P6 because pure states are precisely those states, p for which 
p 2 = p, and this fits in nicely with P6. 

In order to state a precise result, let us introduce the tensor dku obtained from 
considering the Jordan structure of the Aj. Specifically, if {X k , A;} denotes the Jordan 
commutator of Afc, A/, then 

4 d2 ~ l 
{Afc, A;} = AfcAz + AjAfc = -Jd^M + dkn\i 



We use the dku to introduce an operation amongst vectors x, y G R d2 ~ l , via 

d 2 -l d 2 -l n 2 -l 

x U y = ( ^2 dijkXjVk, d 2j kX j y k , . . . , ^ d ijk Xjy k , . . .) 



x U y is thus a vector in R . We can now state 

Proposition 4. Every density matrix can be represented in the form in Equation 
(EH) with (3 = 2s/3b + ^u*, where & is any vector in Bf' 1 with || #, || 2 < f and 

K = +y d2 ~ 2 ^ ^ ■ Conversely any Hermitian matrix admitting such a representation 
is necessarily a density matrix, p is pure precisely if it can be represented in the form 
in Equation EU with < (3, (3 >= and (d - 2)/3 = (3 U (3. 

The proof is straightforward. Since p = H 2 and H itself can be expanded as a linear 
combination of Id and the Aj (albeit with the coefficient of Id different from ~), the 
first part of the result follows from the linear independence of {Id, Aj}. For the second 
part, we represent p as in Equation (|3.1|) and equate it to its square. 

Once again, the difficulty is in the analysis of states which are not pure. It is worth 
mentioning that the pure state condition in the above proposition is essentially the 
same as that obtained from the characterization P5 (for a pure state the characteristic 
polynomial is p(t) = t d + (-l)^ -1 , i.e., b x = 1, 6, = 0, % > 2). 

It should be pointed out that even an analysis of the pure state conditions is far from 
trivial. The condition (d — 2)j3 — j3 U (3 is vacuously true when d = 2 (since the Pauli 
matrices anti-commute) . For d > 3, this condition imposes genuine restrictions. It is 
an interesting problem to find an orthogonal basis for the space of Hermitian matrices 
(the generalized Gell-Mann matrices form just one amongst many) which is close to 
"abelian", i.e, one for which many of the d^u vanish, to facilitate the analysis of the 
condition (d - 2)(3 = (3 U (3. 

In contrast, the parametrization discussed in the next section yields a very simple 
characterization of pure states. 



4. A DIFFERENT PARAMETRIZATION OF POSITIVE MATRICES 



In this section we recall informally the main result of 5J on the parametrization of 
positive matrices. In order to do that a few preliminary definitions and notions are 
needed. 

To any contraction T, one defines its defect operator via 

D T = (I - T*T) 1/2 

Here M* is the adjoint of an operator (when M is a scalar, this is merely complex 
conjugation). 

To such a contraction one can also associate a certain unitary operator, called the 
Julia operator of T via 



(4.1) U(T) - 

Thus, U (T) is a unitary dilation of T. 



T D T * 
D T —T* 



If we are given a family of contractions T k j,j > k with Tkk = for all k, then we 
associate to it a family of unitary operations via the Julia operator construction as 
follows. We first let U k>k = Wj while for j > k we set 

Ukj = Uj-k(Tk t k+i)Uj-k(Tk t k+2) ■ ■ ■ Uj-k(Tk,j)(Uk+i,j © Iv r * )j 

where 

Uj- k (r ktk+l ) = i ®u{r k>k+l ) © /. 

To a family of contractions, T k j one can associate a row contraction via 

flkj = [ r k,k+i, D T * k+ r k>k+2 , . . . , £>r* fc+1 • • • D ri ._V kJ ] 

and a column contraction via 

= [ Tj-ij, 1", 2jA\ j ;,- T kij D rk+1J . . . D Tj _^ ] , 

where " f ' stands for matrix transpose. For more details on the ranges and domains of 
these operators see 0. 

Then the main theorem regarding positive matrices can be stated informally as 
follows (for a precise statement, especially concerning the ranges and domains of all 
operators involved, see [5]) 

Theorem 5. The matrix S = [S k) j] k - =1 as above, satisfying S* k = S k j, is positive if and 
only if i) S k k > 0, k = 1, . . . , d and ii) there exists a family {T k j \ k, j = 1, . . . , d, k < j} 
of contractions such that T kjk = for k — 1, . . . , d valid, and 

(4.2) Sk,j = L* kjk (R k>j _ 1 U k+li j^ 1 C k+ljj + D T * k k+i . . . D T * ._r ktj D rk+l j . . . I) v . 

where L kjk is any square root of S kk . 

Definition 6. The contractions T k j, with j > k, will be called the Schur-Constantinescu 
parameters of S. 



These parameters were first discovered for Toeplitz matrices by Schur, [23], albeit 
in the guise of a problem about power series which are bounded in the unit circle. 
In our humble opinion, it was our late colleague and teacher, T. Constantinescu, who 
championed the study of these parameters to cover all positive matrices (more generally 
to matrices with displacement structure, jS]) and most adroitly brought to fore many 
of their interesting features. Therefore, we have chosen to call these parameters, the 
Schur-Constantinescu parameters, in his honour. 

We will illustrate Theorem (JHJ) via the case of 3 x 3 positive matrices. 

Oil D12 Sl3 

Thus, let S = 5*2 522 be a positive matrix. Then Sa > and let us pick 

. ^13 ^23 ^33 _ 

La as the positive square roots of Sa. In this case L* u = La. Then per Theorem 
there are complex numbers Ti 2 , Ti 3 , T 2 3 in the unit disc such that gives: 

5l2 = L{ 1 Ti2L22, 

523 = -^22^23-^335 
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S13 — l* u (r 12 r 2 3 + -D r * 2 r 13 .D r23 ) l 33 , 



Note that there is a recursive procedure to determine the T k j. The first and the 
second equations yield ri 2 ,r 2 3 from quantities already known, while the last equation 
yields T 13 from quantities already determined at the first two equations. 

Whilst, the Schur-Constantinescu paramters are defined directly in terms of the 
entries of S, one could also seek expressions for them in terms of the vector (3 of 
Equation (|3.1|) (i.e, when S is a density matrix). See [7J |3] for such expressions. In 
particular, for d = 2 the analogue of the Bloch sphere is now a cylinder. 

It is appropriate to make several comments about these parameters at this point: 

• CI As can be expected from the form of Equation ()4.2|) . Theorem (jSJ) is valid 
for operator matrices, i.e., matrices whose entries are matrices or even operators 
in infinite-dimensional spaces, i.e., for elements of ® C(H), with TC allowed 
to be infinite-dimensional. In fact, one can easily extend the result to infinite 
matrices with (possibly infinite-dimensional) operator entries. 

• C2 Though we only called the as the Schur-Constantinescu parameters, a 
full parametrization is provided by the rf ^ ^ contractions T k j, k < j and the 
L iiy i = 1, . . . ,d. In the case of scalar valued matrices, i.e., when 7i = C, we 
thus get the right count of d 2 real parameters. Note the T kk = are just some 
fake parameters, included in the statement of the theorem to avoid an artificial 
separation of the j — k + 1 case from that for other values of j. 

• C3 Since the La are allowed to be any choice of square root of Su (i.e., 
Su = LuL^), the parametrization will be different for different choices of the 
La. A most natural choice would be the Cholesky factorization of Su. In fact, 
as described in [3] , there is an algorithmic proof of Theorem which automat- 
ically yields the Cholesky factorization of S. In the infinite-dimensional case, 
some of the algorithmic flavour of the proof is lost. 

• C4 While, Equation (J4.2)) in Theorem (jHJ) is nonlinear and looks quite com- 
plicated, there is an iterative feature to it (as mentioned in the 3x3 example 
given before), inasumch as in each equation there is just one of the T k j being 
solved for. It is precisely because of this that the Schur-Constantinescu param- 
eters have an inheritance property, namely that the parameters of any leading 
principal submatrix (recall these will be positive themselves) are the same as 
that obtained from the original matrix. 

• C5 Since the proof of Theorem © supplies the Cholesky factorization of S, 
we get an algorithmic recipe for finding one Kraus operator representation of a 
quantum channel $. Since the Cholesky factor, V is lower triangular, the Kraus 
operators, Vi, thereby obtained from V (as described in Section 2), tend to be 
sparse. This can be useful in determining sufficient conditions for a channel to 
be entanglement breaking, or for computing quantities associated to channels 
such as the entanglement fidelity, for instance. The utility of using the Cholesky 
factorization lies not just in the avoidance of spectral calculations (as would be 
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the case if T was found from the spectral factorization of S), but that most of 
the Kraus operators Vi are then sparse. 

• C6 Returning to a positive matrix, S, whose entries are scalar, it is known that 
if Su = 0, for some i, then the entire row and column to which Su belongs 
has to be zero. Therefore, a reasonable convention to assume is that Tkj = 0, 
whenever SjjSkk = 0. With this convention, the Tkj-, La provide a one-one 
parametrization of positive matrices. 

• C7 In the previous section we saw that even the problem of characterizing pure 
states via the proposed parametrizations of that section was not fully resolved. 
However, the Schur-Constantinescu parametrization provides a very simple and 
effective characterization of rank one states, viz., S is rank one iff all T^j = 0, 
except for those cases in which SjjSkk 7^ 0, in which should be on the 
unit circle. 

• C8 Let S be positive matrix. Then there is a very simple formula for its 
determinant in terms of the Tkj, viz., 



This is useful since some entropic quantities can often be expressed in terms of 
determinants, [T5] . 

• C9 While Equation (|4.2|) is intricate, there is a useful diagram (called a trans- 
mission line diagram) which keeps track of all the matrix products in it. 



In this section two additional applications of the parametrization of the previous 
section are provided. The first is to show that block Toeplitz states have positive 
partial transpose. The second is to examine the restrictions on the relaxation rates for 
an open quantum iV-level system imposed by the requirements of complete positivity 



5.1. Toeplitz States. The positive partial trace condition of |2lH IT2*] has been found 
to be a very useful operational condition for entanglement. While, for general states, it 
is known to be necessary and sufficient only for 2x2 and 2x3 states, there have been 
several arguments in favour of the notion that states which satisfy this positive partial 
trace condition (PPT states) are "close" to being unentangled, at least inasmuch as 
they are not useful for tasks such as dense coding. Similarly there have been several 
attempts at studying the PPT property for positive matrices which satisfy additional 
conditions, see[2j. In this section we provide a contribution along the same vein. We 
show that positive Toeplitz matrices are PPT states. 

The proof of this result was first found by considering the Schur-Constantinescu 
parameters for 3x3 block Topelitz psoitve matrices. This proof can be extended in 
a simple but tedious manner for x d 2 states. But there is, in fact, a second proof 




5. Two Further Applications 



(cf, usa). 
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which works for all dimensions. We provide this first and then discuss the parameter 
based proof. 

Proposition 7. A Toeplitz mixed state is PPT. 
Let A e C NxN be a Toeplitz matrix given by 



a a_i 
ai a 



We will first, for illustration purposes, show that A T is also positive. This is, of 
course, true for arbitrary positive matrices, but it will serve to illustrate the proof in 
the partial transpose case. Then the ij-ih entry of A is given by A^ = cii-j. The 
transpose of A, denoted by A T , is 



a ai 
a_i a 



' ' 0"a-l 



a -n a -n+l ' ' ' 

,with Ajj = a,j-i. Next, in the cycle notation, let a"o be the element of the symmetric 
group Sn on iV letters, {1, 2, • • • , iV}, defined by 

a = H (k(N-k)). 

l<k<N 

cr induces two simple operations on iV x N matrices. If M e C NxN takes the form 

Wi 



M 



U>2 



W, 



, where w^s are rows of M, we define the operation R ao by 



, i.e. R ao simply permutes the rows of M as specified by oo- Another operation on 
columns, C ao , is define in the same way. Now we notice that if A is Toeplitz as given 
above, then 



[R ao (C ao (A))} id - Ajv-i,jv-j - aj-i - Afj. 
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Since R ao and C ao preserve the characteristic polynomial, we have shown that a if a 
Hermitian Toeplitz matrix is positive then so is its transpose. 



The above fact can be extended to the partial transpose of an NM x NM Toeplitz 
matrix A in the following way: Let a m be the same permutation as a"o on the letters 
{run, ran + 1, • • • , (m + l)n — 1}. If a G S^2 is defined to be the disjoint product 
o"oO"i • • -<Jm-i, and R a and C a are the induced operators, then by the same argument 
as above, we have R a (C a (A)) = A PT , where A PT denotes the partial transpose of A. 
Once again these operations preserve the characteristic polynomial for Toeplitz matri- 
ces and hence if A is positive, in addition, we find that so is A PT . Thus a positive 
Toeplitz matrix is PPT. 



The Schur parametrization of positive matrices gives another proof of proposition 
1 that is immediate. If B is a block Toeplitz matrix, then B is also Toeplitz. So let 
B be, for instance, a 3 x 3 block Toeplitz matrix. Using the Schur-Constantinescu 
parameters and the block Toeplitz property of B, we can write B explicitly as 

Ah(T\ + D r *T 2 D Vl )A^ 
Ahr^ 



A2TIA1 



AzTxAi 
A 



_ A^((T\f + D Tl V* 2 D Tl )Ai AiY\Ai 

, where each entry is an N x iV matrix. Note that due to the block- Toeplitz nature 
of B its Schur-Constantinescu parameters need be indexed by only one subscript. 
Transpose block- wise gives us A PT . By the spectral theorem, D t t = (DfJ T . So, sim- 
ply by inspection, we see that A PT has Schur parameters {(A^) T , Tj , Tf}. Therefore 
A PT > 0. This is in fact true in general: 

Proposition 8. If A e C MNxMN is block Toeplitz, then A is PPT. 

The basic idea is to show that If A is parametrized by {Tj}, then A PT is parametrized 
by {rf}. Note that the block- Toeplitz property means that Schur-Constantinescu 
parameters of A depend only on one index (cf., the 3x3 block case). We will omit 
the proof, which is straightforward but tedious. Via the combinatorial structure of the 
Schur parameters, one can see how the parametrization of A gives rise to that of A PT . 
The so-called " lattice structure" of the Schur parameters for the 4x4 case is shown in 
Figure 1 below. Each transfer box in Figure 1 describes the action of the Julia operator 
U(Ti). 

Let U T (T) denote the transpose of the Julia operator of T, i.e. 

r T (D r f 
(D r »f -(r*f 



[D r r) 



D-pT* 

-(T T y 



U(T T ). 



Each entry of the positve semidefinite kernel {Aij} corresponds to those paths in the 
diagram that startfrom Ljj and end at La*. For example, each path from L 33 to L\\ 
describes to a summand in the expression for A 13 . So we can see that the transmission 
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line diagram of A PT is then obtained by replacing each U(Ti) transfer box by that of 
U T (T t ). 




Figure 1. Lattice structure for 4x4 positive matrices 



5.2. Constraints on Relaxation Rates. In this subsection we revisit the very in- 
teresting work of [22] on the constraints imposed on the relaxation rates of an open 
iV-level quantum system by the requirement that its evolution be completely positive. 
In order to keep the notation the same as in [22]? we will, in this subsection only denote 
the Schur- Constantinescu parameters by gij (and not T^). 

Let us first briefly review the contents of [22]- Let p(t) be the state of an open 
iV-level quantum system and let p be the vector in C N ' 2 which represents vec(p). Then 
its evolution can be expressed via 

(5.1) p={- l -L H + L D )~p 

n 

where Lh and Lo are N 2 x N 2 matrices representing the Hamiltonian and dissipative 
parts respectively of the evolution of p. Let the index (m, n) denote the number 
m + (n — 1)N. Then the non-zero entries of Lp are given by 

\Ld) (rn,n){rn,ri) L mn , 771, ^ Tl 

(LD)(m,m)(l,l) = Imhm^l 

N 

(-^d) (m,m)(m,m) ^ ^ Ikm 

Here 7& n is the population relaxation rate from level | n > to | k >. The jkn are real 
and non-negative. Tk n (for k ^ n) is the dephasing rate for the transition from | k > 
to | n >. Since, T kn = T nk , it is easily seen that (L D )( m>n )( mjn ) = {L D )(n, m )(ri,m)- A key 
step in the work of [22J is to express Tk n as a sum of two summands, in recognition of 
the fact that dephasing is also enhanced by population relaxation, to wit 

v -r p j-H 

1 kn — 1 kn' 1 kn 

with T p kn , the decoherence rate due to population relaxation and Tf n the decoherence 
rate due to pure phase relaxation. The requirement that the open quantum system's 
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evolution be completely positive, [TJ Ej, imposes restrictions on 7^ n and Tk n - These 
restrictions can be expressed as the requirement that a certain (iV 2 — 1) x (iV 2 — 1) 
matrix concocted out of the 7^ n and Y kn be positive, [22] ■ However, per [221 , this 
requirement can be reduced to verifying that a related (N — 1) x (N — 1) matrix be 
positive. The form of this (N — 1) x (N — 1) matrix will depend on a choice of an 
orthogonal basis for the space of traceless, Hermitian (iV 2 — 1) x (TV 2 — 1) matrices. 
However, positivity of this matrix itself is independent of the choice of basis. The 
excellent analysis of [22] is unfortunately marred for the N = 4 case by an incorrect 
criterion for positivity. Indeed, Equation (28) of [22] are only necessary for positivity, 
while Equations (31)-(32) are (as correctly claimed in (221) also just necessary (though 
they come closer to sufficiency than Equation (28) of [22])- 

In the sequel, we will use the Schur-Constantinescu parameters to analyse the N = 4 
case of [22] • As in [22] the evolution Equation (j5.1j) is completely positive iff the 3x3 
real symmetric matrix B = (by) is positive. To specify the entries of the bij, we denote 
by Tf ot the quantity \ Y^n=2 El^i ^mn- Then the entries of B are given by 



foil 


pd 
~ 1 tot ' 


_ (yd , yd 
\ L 13 ^ 1 24 




-pd 
~ 1 tot ' 


_ /pd — 1 — yd 
V 1 13 > 1 24 




-pd 
~ 1 tot ' 


_ (r d _|_ W 

I 1 12 ^ 1 34 




(r? 2 


1 34J 

2 


&13 


(it 


-rf 3 ) 

2 


fin, 


(r? 3 


1 24J 



Now i? is positive iff bu > 0, i — 1, . . . , 3 and the Schur-Constantinescu parameters 
9i2,gi3,g23 are in the closed unit disc. Since B is real this is equivalent to demanding 
that the g^ belong to the interval [—1, 1]. 

The conditions bu > become 

-pd _j_ -pd _j_ -pd _j_ pd \ -pd _j_ pd 

12 14 23 34 — 13 24 

pd _i_ "P^ J- J- > p^ _i_ p<^ 

12 13 24 34 — 14 23 

pd I pd I pd I pd \ pd I pd 

13 14 ' 23 24 — 12 34 

Now 6i2 = v^nfi f i2V / ^22- So g u G [—1,1] becomes 

4r^r^ 4 -(r?3-rfj 2 -(r^3-r^3) 2 +(rf3-r^ 4 ) 2 +(rf 4 -r^3) 2 -(r? 4 -r^ 4 ) 2 -(r^3-r^ 4 ) 2 >o 

Likewise the condition #23 £ [— 1, 1] becomes 

4r? 3 r* 4 -(r? 2 -r? 4 ) 2 -^ 

Finally gi 3 G [—1,1] becomes 

bllb 2 2b 3 3 + 2&12^13^23 > ^11^23 + 6 22^ 3 + &33&12 
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Note that the condition | |< 1 is not similar to the condition for the other g^j to 
be in [—1, 1]. This is to be expected since the formula for g^ for k > j '• + 1 is more 
intricate than those for the gjk, k — j + 1. Furthermore, this last condition is precisely 
one of those obtained in [22J. However, the conditions obtained here are necessary and 
sufficient. 



6. Conclusions 

Since positive matrices play a vital role in many applications, it is of importance 
to obtain computable parametrizations of them. In this paper we discussed several 
such potential parametrizations. Which one of them one ought to use is, of course, 
a matter dictated by the application one has in mind. We argued, hopefully persua- 
sively, in favour of the versatility of the parametrization proposed in There are 
several other applications besides the ones discussed here, to which one could apply 
this parametrization. This will be the subject of future work. 
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